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Abstract 

The spatial gradient expansion of the generating functional was recently devel- 
oped by Parry, Salopek, and Stewart to solve the Hamiltonian constraint in Einstein- 
Hamilton- Jacobi theory for gravitationally interacting dust and scalar fields. This 
expansion is used here to derive an order-by-order solution of the Hamiltonian con- 
straint for gravitationally interacting electromagnetic and scalar fields. A conformal 
transformation and functional integral are used to derive the generating functional up 
to the terms fourth order in spatial gradients. The perturbations of a fiat Friedmann- 
Robertson- Walker cosmology with a scalar field, up to second order in spatial gradients, 
are given. The application of this formalism is demonstrated in the specific example 
of an exponential potential. 



1. Introduction 

Hamilton- Jacobi (HJ) theory has many apphcations in the perturbative and non-perturbative 

analysis of dynamical systems in classical mechanics. Peres developed an Einstein-Hamilton- 
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Jacobi (EHJ) formulation of general relativity in which a generating functional has to satisfy 
the momentum and Hamiltonian constraints of general relativity In the framework of 
quantum cosmology, it was known that the momentum constraints require any wave func- 
tional to be diffeomorphism invariant 0. In a WKB approximation, such a requirement 
translates into the diffeomorphism invariance of the generating functional. The Hamilto- 
nian constraint is a non-linear functional partial differential equation that governs the time 
evolution of the generating functional. 

Based on Peres' formalism. Parry, Salopek, and Stewart used a series expansion of the 
generating functional in spatial gradients of the fields to derive an order-by-order approx- 
imate solution of the Hamiltonian constraint for general relativity with matter fields (see 
Ref. m, from now on referred to as PSS). Such a generating functional is diffeomorphism 
invariant in each order of the expansion. Salopek and Bond used this formalism to show how 
non-linear effects of the metric and scalar fields may be included in stochastic inflationary 
models. The main advantage of this analysis is that the lapse function and shift vectors do 
not appear in the EHJ equations. Therefore, one obtains a coordinate-free approach to cos- 
mological perturbations. In the above models matter fields consist of self-interacting scalar 
and dust fields. 

In this article the above formalism is extended to minimally coupled gravitationally 
interacting scalar and electromagnetic fields. Such minimally coupled electromagnetic fields 
give rise to conformally invariant field equations. Hence, the electromagnetic field energy 
density is proportional to where a is the scale factor. Consequently, the electromagnetic 
field is diluted away during the De Sitter expansion phase of the inflationary cosmologies. 
To break the conformal invariance, a direct coupling of gravity to electromagnetism or 
corrections due to the quantum conformal anomaly have been considered 0. A coordinate- 
free approach to the perturbative analysis of cosmological models with electromagnetic fields 
could eventually lead to a better understanding of the primordial magnetic fields. The 
generating functional up to the third order in the spatial gradient expansion is given in 
section ^. Following PSS, section ^ is a demonstration of how a recursion relation and a 
functional integral in superspace can be used to derive the higher order terms in the spatial 
gradient expansion from the previous terms. Section ^ is an exhibition of the gauge fixing 
and the solution of the field equations. 
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2. ADM reduction and the EHJ equations 



The action for minimally coupled gravitationally interacting neutral scalar and electromag- 
netic fields can be written as 

C /(4) (4) 

X = J V 9[\R -kM'^ - V{(P) - \F,,F^^]d'x, /X = (0, 1, 2, 3), (1) 

where F^^, = A^^^ — A^ ^ is the electromagnetic field strength and Vi^cj)) is the scalar field 
potential. ADM reduction of the above action is achieved by defining the 3-metric 7^1, := 
Qnu + n^rii, and the vector potential A = A|| + Aj_, such that (A||)^ = 'j^Ai,, (A^)^ = 
— {n'^ A^)nfj^, where is the unit vector field normal to space-like hypersurfaces of simul- 
taneity parametrized by t. In the basis {dt,di),i = 1,2,3, such that n^{di)^ = and 
n^' = {dt, -N'di)/N (no sum), the following relations hold: (A||)i = 7^^^, Aq = NA^+N'Ai, 
and 

ds^ = -N'^df + -f^j{dx' + N'dt){dx^ + NUt). (2) 

Then one follows with the procedure outlined in Ref. to derive the Lagrangian L for gravi- 

tationally interacting electromagnetic and scalar fields. With the momenta vr*-' = , Tr*^ = 

SL 8 L 

— E'^ = — where ' := d/dt^ after a Legendre transformation the Hamiltonian is 



5A, 



Hamiltonian = J {Ni^Hf, + Aog)d^x, (3) 



where 



Ho = 7-'/'vr^%'='(27,,7,fc - 7,,7h) + 7'/'[^(0) - i?/2 + F^'Fu/4 + ^(^^2] 

+ 7-1/2 [^i^. _^ {■n't>f]/2 = 0, Hamiltonian constraint. 
Hi = —2nly + FiiE^ + n'^4>A = 0, momentum constraint. 



(4) 



-E'^\i = 0, Gauss law constraint. 



|j is the 3-covariant derivative (for covariant derivatives of tensor densities and sign conven- 
tions see PI). Utilizing Hamilton's equations, the evolution equations for the fields are as 
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follows: 



^ = N-f-^/\^ + N'(f),i, (5) 

= N2-f-'/\''\2-fiajk-1ij1ki) + 2N^i\j), (6) 

A = N^''/'E, + N^Fj, + Ao,.. (7) 

The evolution equations for the momenta are considerably more complicated. They are given 
by 



TT 



- F\P^ + If^FimF'"^ + i7*'0,^0'' - rr + 7'''^(0) } 

+ ^7^/2 (a^I^^' - Ar|/I'7^^') + (A^^vr*^') ^ - 27r'^(^Ar^'\„, (9) 

= 7^/2 i^N^rnF""' + iVF™|^) + (A^™^' - A^'^") ^ • (10) 

Instead of solving the evolution equations for the fields and momenta, one can try to 
solve the EHJ equations. The EHJ equations are derived by the substitutions 

in TC^ and Q. S = Sl'jij, (p, Ai] is the generating functional (Hamilton's principal function) 
0. The Hamiltonian constraint is a hyperbolic functional partial differential equation for S. 
After solving the EHJ equations, (^0) and (|n|) yield the full set of the evolution equations. 

3. The spatial gradient expansion and the order-by-order solution 
of the EHJ equations 

The momentum constraint implies that the generating functional is diffeomorphism invariant 
One such-diffeomorphism invariant quantity is S* = / 7jj, Aj]7^/^(i^x. More gen- 
erally, a diffeomorphism invariant S can be a multiple integral of some multi-point functions. 
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The contribution of such highly non-local terms could be important, for example, if the spa- 
tial inhomogeneities are correlated. However, at least in the lowest orders, the contribution 
of such terms to the generating functional are expected to be insignificant in the generic case. 
Likewise, the Gauss law constraint implies that S is gauge- invariant, e.g. 5* = >S'[-Fjj]. Other 
gauge-invariant quantities like ^ Aidx^ could also be included in 5*. However, if the space-hke 
surfaces of simultaneity are simply connected, one can write all such quantities in terms of 
Fij using Stokes theorem. Non-simply connected three-manifolds are not considered here. 
The Hamiltonian constraint determines the time evolution of the fields. Following PSS, an 
order-by-order solution of the Hamiltonian constraint is achieved by the expansion of the 
generating functional in spatial gradients: 

oo 

S = ^A"5("\ (12) 

n=0 

where A" denotes the number of spatial gradients in S^'^\ It turns out that the scalar field 
is indispensable in this model and dominates the dynamics of the space-time. Moreover, 
in the limit i? ^ 0, in the absence of electromagnetism, the zeroth order solution is exact. 

Therefore, merely based on dimensional grounds, this heuristically suggests the expansion 

R - 
parameter should obey A oc y^^-^ ■ -^^^ ^^^^ relation, R is an appropriate combination of 

the curvature invariants of dimension such that in the flat space limit where the three- 
curvature is vanishing ^ = 0, and V{(j)) represents the potential energy density of the scalar 
field. The convergence of the above series is an unsolved problem [pT|] . 

An order-by-order solution of the EHJ equation is achieved by substituting ([T^) in 
the first equation in (^) and the subsequent expansion of TCq in spatial gradients: 

oo 

Ho = ^A"7i:(") (13) 

n=0 

and requiring the EHJ equation to vanish at each order. In the above equation 

n^'' = 7-^/^^f^(2w. - + 7^/V(0) + 7-^/^(^)72 = (14) 



One can easily obtain the first few terms in ([12D by an ansatz. The zeroth order term 

= -2 y j'/''H{cj>)d\ (15) 
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called the long- wavelength approximation (LWA), is the same as in PSS for some arbitrary 
function H((f)). Inserting S^^^ in ([1^) yields 

-3if2 + ^(0) + 2(^)2 = 0. (16) 

Electromagnetism has no dynamical degrees of freedom at this order. The LWA is very 
important in structure formation after inflation. Therefore, it is unlikely that electromagnetic 
fields play a significant role in structure formation in this model. For n > 

z/W = 0, for n ^ 2. 

(17) 

At each order, = is a linear hyperbolic functional differential equation in the unknown 
functional S^'^\ 

S^"^^ is the next non-vanishing term given by 

S^'^ = J l''\J{m + ^(0)0.0" + L{(t>)F,,F'^)d?x. (18) 

The above terms are not the only terms quadratic in spatial gradients in S^'^\ However, 
the remaining terms are either equal to the above terms modulo surface integrals or vanish 
identically. For example, / '-^^/'^F'^^ eijk(t)'^d?x = — J -^^^"^ {F^^^^eijkcfyd^x = due to Maxwell 
equations -^[jji/c] = 0. It is easy to verify that S^"^^ satisfies the momentum and Gauss law 

5^(2) 

constraints. One notices that ^--7— is already quadratic in spatial derivatives and does 

oAi 

not appear in the second order EHJ equation 7^^^^ = 0. After inserting S''-^^ in the second 
order EHJ equation and grouping together the coefficients of R, 0''0,«; and Fi^F^^ one 
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respectively has 

dHdJ 1 <fj dHdK 1 

dJ dH dHdL 1 

At first sight, (|19|) seems to be an over-determined system for three unknown functions, J, 
K, and L. However, by solving S2 = for K, one can show that S3 is not independent and 

ds2 dsi f dH\~^ 

obeys the relation S3 = 2— — si H — —H —— . In the spatial gradient expansion of 

d(p d(p \d(j) J 

the generating functional 5" for the scalar fields in the absence of electromagnetism, there is 
no contribution from odd order terms. Electromagnetism makes non-trivial contributions to 
odd orders. The only non-vanishing term in S^^^ is 

S^^'^ = j j'/^M{(P)r^ye,kiF'''d^x. (20) 

All other third order terms like F^^\j(j)i, F^^F^^emni'P.p) F^^\ij either vanish or are total 
divergences. To solve the third-order EHJ equation Ti^^-* = one has to use the relation 
r*^ = 7~"^^^(7"'^^^),i. The solution yields 

dMdH 

MH + = 0. (21 

non-trivially satisfies the momentum and Gauss law constraints. To show this, one 

frequently uses the identity Vseijk = Viesjk+Vjeisk+Vkeys for any vector field Vi, achieved from 
S 6 

^gijfc^ = — (Yh^^l'^^'^imn) and multiplication of both sides by Vr- Equations ([T6|) , 



(pj|), and (^1]) form a set of differential equations, easily solvable for most relevant potentials. 
Nevertheless, the full set of differential equations becomes increasingly complicated at higher 
orders. As in PSS one can use the expression for 5**^°^ and the conformal transformation 

= 7.,fi-'(«), u := j (-2^)-M0, ^ = Hn, (22) 

to solve the EHJ equations (^^. The EHJ equations transforms into 

5uix) 



I fij 



n^""^ (23) 
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where 



p=i 

~(2) ^ /V2{^/2 + A.^Hn)uiU,P - l^F.kFimff"'}, z>(") = 0, for n ^ 2. 
R is the conformal curvature and all indices are raised and lowered with the conformal metric 

fij- 

The proof of the integrability of (ESp for gravitationally interacting electromagnetic 
and scalar fields proceeds similar to what was outlined in ||T^. By using the expression for 
^(2) from (0), a functional integration of (|23| ) gives rise to the following expression for S**^^^: 

^(2) = j /1/2 (^j(u)R + k{u)u^^u'' + l{u)FikF'^'^ £x, (25) 



where 



The complementary functionals D' are constants of integration. In the next order — - — = 

ou 

0, therefore S^^^ = J f^^'^F^^ ■jeikiF^''d?x is the most general form of S''-^-* in which ; and Sijk are 
the covariant derivative and Levi-Civita tensor associated with fij respectively. Conformal 
transformation of this expression gives rise to 

5(3) = j Y^n^r=\^e,uiF^U^x. (27) 

Functional integration of p^ ) in the next order gives rise to 



(28) 



5(4) = J d^xf^/^{-i{u)R'^Rij - (3£(m)/8 + m{u))R^ - n{u){R''J - pR/2)u^^Uj 

+ r{u)u'lu'^Uj + s{u){FiraF^'^f + t{u)u^pU'PFimF^"' + v{u)u^rnU'''F"''Fni 

+ w{u)u,^F^'F'\,n + x(m)F-V„F-, + y{u)FkmFrF^^FK^ + z{u)RFinF'^ 

+ a{u)RuF>^^Pm}. 
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where i{u),m{u),n{u),r{u) are given in PSS. In the above expression s{u),--- ,a{u) are 
defined as 

J TT 

s'{u) = Q-^[-llP/A + {—)-\Q-^/A + a\u) = 8jlQ-^ 

z'{u) = (^-^jl + i(^)-2(l/8 + Hm/2 - Hjn-^/A - H^jl)^ , x'{u) = -SPQ-\ 

v'{u) = Q-\-AlH + 3^-72), y'{u) = 8l^Q-\ 

w'{u) = Am-^{21H + 

(29) 

in which ' := d/du. 

4. The evolution equations of the fields 

Once the EHJ equations are solved, the evolution equations for the fields are obtained from 
and (|ll]). A judicious choice of gauge greatly simplifies the field equations. In the 
almost synchronous gauge (A^* = 0) if m is the time parameter, from it follows that in 
the LWA (i.e S = S'^'^^) the lapse obeys A^'-"'^^ = 1. The superscript (n) means that the right 
hand side of the equation contains terms up to {n — l)th order in spatial gradients. 

The choice of u as the time parameter is valid as long as the geometry is sufficiently 



close to that of the homogeneous models (for a relevant discussion see Ref. |]T3[)- Then it is 
useful to replace (H) with the equivalent evolution equation 

4 = 2Nn-'f-'/'-^{2f^Ju - f^,fkl) - 2Hf,, (30) 
for the conformal metric fij which is related to jij via 

7,, =exp|- j{^r'Hdc^]k- (31) 

In the LWA /^^^^■* = 0. f^J^ is the seed metric that contains no dynamical degrees of freedom. 
The first non-trivial evolution equation for Ai in the temporal gauge Aq = 0, is obtained 
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from the conformal transformation of (|^: 



A] 



(3) 



6S^ 



1 17m 



(32) 



Since d/du and ; commute, the evolution equation for Fij is easily derived from the above 
equations to be 



[j;i]m-: 



with the solution 



F^f = exp |- (^8 y du n-H{u)^ 5g55]VfcV"^| 



(33) 



(34) 



In the above equation and what follows, J^ij is an arbitrary antisymmetric tensor field, V 
refers to the covariant derivative with respect to the seed metric and the indices are raised 
with the seed metric. The exponential of the matrix differential operator is defined as: 



exp {■ ■ ■ } := |(5r55 - y Q-H{u) du^ S^/^^VkV"" (35) 

+ ^ (s y n-H{u) du^ VpV^^f/^] VfcV'" + ■ ■ ■ I ^^z. (36) 

Once the evolution equations for the fields are solved, the evolution equations for the mo- 
menta are easily derived from ([TT|). In particular, the electric field obeys the equation 



6Ai 



-Af'\u)F'^',, 



In higher orders (^) shows that A^(") ^ 1. For example. 



du 



dl 

du 
(1) 



dH\ 



(37) 



(3^ 



where B!^^^ is the three-curvature associated with the seed metric fl^ . Obviously, the higher 
order evolution equations are more complicated. The third order evolution of the conformal 
metric is derived from (|30|) and ( ^8|) to obey 



p(3) 



I TiiT f{l)kmf{l)lnf{l) 
~ ^ kl^ ran J J J ij 



H dH 2f TT d.j ^ 



H ^dH ^ '^^ ^ _j_ 2^ 

2 d(h du 



- S/^i^^.^/^i)™" 1^ . (39) 
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As a demonstration of an application of the formalism developed so far, one could 
compute flj^ for an arbitrary magnetic field J-'f^i and seed metric f^j^ with a scalar potential 



V = Voexp l^y'^'^j'' "^^^ general solution of (|T6|) for p 7^ 1/3 is given in ||T^. The general 
parametric solution of ([T6|) with H and as functions of an independent variable v is, 



H 



Vo 



2 /3 
exp{— -} coshf , 



+ 



±2 + ^ 



(40) 
(41) 



where 0m is the integration constant. A special solution of (|T6|) for p 7^ 1/3 is 

1 



H 



exp{^=}, 



that corresponds to the Halliwell attractor [Q. By using (^), (|5UD, ( |1^ ) and with a choice 
of time parameter such that lim^^o^ = —00, the contribution of all the terms second order 
in spatial gradients to the evolution of the conformal metric at this order is 



/(3) 
Jij 



-2c- 



+ 



(p + l) [piSp-1] 



u 



-2p+l^(l) 



P) 



1 -2p 



p(3p 



u 



[-\:FmnJ''^^flf + 2^,"^,„) , p ^ 1 (43) 



where c is the integration constant and the seed metric is absorbed to raise indices. Integra- 
tion of the above equation yields 



(3) 



p)(l — 2p) \_p{?>p 



-2p 



U 



-4p+2 / 1 'pmn _i_ -77 " -77 \ 

1 4«' mn'' J ij ~T Ji J jn j 



+ 



P^ — 1 

After taking (|3TD and (^ 



(3) 
lij 



p{3p 



u 



(1) 



(44) 



into account, the three-metric is given by 



Vn 



(l-p)(l-2p) 



p{3p - 1^ 
Vo 



u 



-2p+2 ( 1 -77 ^mn r{i) _i_ -77 " -77 



p{3p 



(45) 
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Two points regarding the above equation are in order. First, as a perturbation of a flat 
Friedmann-Robertson- Walker (FRW) model, the growth of inhomogeneities from terms sec- 
ond order in spatial gradients is encoded in ^ylj . One notices that with the choice of u as 
the time parameter, such a growth resulting from the curvature is insensitive to the exact 
form of the exponential potential. In other words, deviations from a flat FRW cosmology 
resulting from a magnetic field evolve quite differently from that of the inhomogeneities due 
to the space-time geometry. In particular, one notices that if p > 1, the contribution of a 
primordial magnetic field to spatial inhomogeneities decays rapidly. 

Secondly, the reader is cautioned against a mini-superspace perturbation of a flat 
FRW model, namely, a perturbation within a mini-superspace cosmology. For example, for 
a vanishing magnetic field in a Bianchi I cosmology where the conformal three-curvature 
is vanishing {i.e. Rij = 0), one would get the obviously incorrect result 'jlj\t) = lifit) 
for n ^ 1 (because S**^"-* = for n ^ 0). The exact evolution equation for such a model 
could be derived from (@-10) and which could be solved numerically for generic initial 



conditions. This problem also exists in the dust model treated in PSS which has been quite 
successful in deriving the contribution of short wavelength fields to structure formation and 
the anisotropics of microwave background radiation. 



5. Conclusion 

The spatial gradient expansion of the generating functional was developed by PSS to solve the 
Hamiltonian constraint in EHJ formulation of general relativity for gravitationally interacting 
dust and scalar fields. The spatial gradient expansion could be consistently apphed to solve 
the Hamiltonian constraint for gravitationally interacting electromagnetic and scalar fields. 
At each order, the EHJ equation is a linear functional partial differential equation in the 
unknown functional S''-"^ which after a conformal transformation could be integrated to yield 
g{n)^ g(2) ^(3) y^QYe calculated in detail and S^^^ was given. Such an order-by-order 
solution of the EHJ equation gives rise to order-by-order corrections to the fields evolving 
in a flat FRW model. The corrections are due the presence of spatial inhomogeneities and 
magnetic field. Not surprisingly, such corrections start with terms second order in spatial 
gradients. The formalism was applied to the specific example of a scalar field with potential 
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f2 

V = Voexpi — 4 Contributions of all the terms second order in spatial gradients to the 

V V 

metric were derived. 

I would like to thank Ali Mostafazadeh, Don Page and Dave Salopek for helpful 
discussions. 
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